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Abstract 

The evolution of pairwise quantum correlations of Bell cat-states under amplitude damping 
is examined using the concept of quantum discord which goes beyond entanglement. A closed 
expression of the quantum discord is explicitly derived. We used of the Koashi- Winter relation. A 
relation which facilitates the optimization process of the conditional entropy. We also discuss the 
temporal evolution of bipartite quantum correlations under a dephasing channel and compare the 
behaviors of quantum discord and entanglement whose properties are characterized through the 
concurrence. 
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1 Introduction 



Quantum entanglement plays an important role in the area of quantum information such as quantum 
teleportation [I], superdense coding [2], quantum key distribution [3], telecloning [5j and many more. 
It is one of the most fundamental concepts in the description of quantum correlations in multipartite 
quantum systems and makes possible tasks in quantum information which are impossible without 
it. Therefore it was important to quantify the amount of quantum correlations in a given bipartite 
quantum state. Much effort has been devoted to the classification of quantum states into separable 
and entangled states (see O [7] and references therein). Until some time ago, entanglement was 
usually regarded as the only kind of nonclassical correlation in a composed state. However, there 
are other nonclassical correlations different from those involved in entanglement which are useful 
for quantum technology. Hence, to characterize all nonclassical correlations present in a multipartite 
system, the so-called quantum discord, which goes beyond entanglement, was introduced in [H[9]. Now, 
it is commonly accepted that quantum discord constitutes a new resource for quantum computation. 
Among the evidences of the relevance of the quantum discord, one may quote for instance quantum 
non-locality without entanglement \10\ \TT\ [T2] and the advantages offered in increasing the rapidity 
of certain computational tasks with separable states in comparison with their classical counterparts 
[13 1 116]. But, it may be noticed that despite increasing evidences for relevance of the quantum 
discord in information processing tasks, there is no straightforward method to get the analytical form 
of quantum discord in a given quantum state. This is mainly due to the fact that its evaluation 
involves an optimization procedure which is in general a hard task to perform. A reliable algorithm to 
evaluate quantum discord for general two-qubit states is still missing and only few analytical results 
were obtained for some particular forms of the so-called two-qubit X-states [TTJ, [181 [13 1201 EH (22] . 
This class includes the maximally entangled Bell states and Werner states [23]. A closed expression 
for the discord of arbitrary states remains an important challenge. 

On the other hand, quantum optical tools are expected to be useful in the context of quantum 
information science, especially for communications using qubits over long distance. However, optical 
qubits suffer from decoherence due to energy loss or photon absorption. To reduce the decoherence 
effects, encoding qubits in multi-photon optical coherent states seems to be a promising alternative. 
They are more robust against small levels of photon absorption (see [23] ) . The photon loss or amplitude 
damping in a noisy environment can be modeled by assuming that some of field energy and information 
is lost after transmission through a beam splitter. 

Considering this problem and motivated by many works devoted to study entanglement properties 
in a composite system involving coherent states (for a recent review see [2S]), we investigate in this 
paper a method to describe the evolution of quantum discord in Bell cat-states under amplitude 
damping. We present an algorithm to calculate the quantum discord. First, we obtain an explicit and 
simplified expression for the conditional entropy, exploiting the Bloch representation of the density 
matrix. Then, we combine the purification method and the Koashi- Winter relation [2U] to perform 
easily the optimization of the conditional entropy. This allows us to get a closed form of quantum 
discord in damped Bell cat-states. Exploiting our algorithm, we examine the dynamical evolution of 
the system under a dephasing channel to compare the temporal behaviors of quantum discord and 
entanglement as quantifiers of non-classical correlations. 
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This paper is organized as follows. Section 2 provides an introduction to quantum discord. We 
present the main definitions and properties. Section 3 concerns Glauber coherent states superposi- 
tions subjected to an amplitude damping channel. In the optical context, amplitude damping can be 
appropriately modeled by having the signal interact with a vacuum mode in a beam splitter. Having 
identified the effect of amplitude damping, we study in section 4 the evolution of the quantum corre- 
lations of Bell cat-states under amplitude damping. We obtain the explicit form of quantum discord. 
Section 5 deals with the comparison of the evolution of quantum discord and entanglement under a 
dephasing channel, using the results of the previous sections. Finally, the last section recalls the main 
results of our work and suggests further issues deserving to be investigated. 



2 Quantum discord: Generalities 

For a state p AB of a bipartite quantum system composed of particle A and particle B, the quantum 
discord is defined as the difference between total correlation I(p AB ) and classical correlation C(p AB ). 
The total correlation is usually quantified by the mutual information / 

I(p AB ) = S{p A ) + S{p B )-S{p AB ), (1) 

where p A ( B ^ = Ttb(A)(p AB ) is the reduced state of A(B), and S(p) is the von Neumann entropy of a 
quantum state p. Suppose that a positive operator valued measure (POVM) measurement is performed 
on particle A. The set of POVM elements is denoted by Ai = {M k } with M k ^ and Mj. = I. 
We remind that the generalized positive operator valued measurement is not required. Indeed, it has 
be shown in [27] that for the optimal measurement for the conditional entropy is ensured by projective 
one. Thus, a projective measurement on the subsystem A project the system into a statistical ensemble 
{p B ,p B }, such that 

n AB . B (M k ®l)p AB {M k ®T) 

p > Pk — b y l ) 

Pk 

where the von Neumann measurement for subsystem A writes as 

M k = UIL k rt: k = 0,1, (3) 

with Ilfc = |/c)(/c| is the projector for subsystem A along the computational base \k), U € SU(2) is a 
unitary operator with unit determinant, and 

p B = tr {M k ®l)p AB {M k ®l) 

The amount of information acquired about particle B is then given by 

S( P B )-Y,Pk S( P f), 

k 

which depends on measurement Ai. This dependence can be removed by doing maximization over all 
the measurements, which gives rise to the definition of classical correlation: 

C{p AB ) = max^ [s{p B ) - Ek Vl S(p B ) 
= S(p B ) - S min 
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where S m i n denotes the minimal value of the conditional entropy 

S = Y,PkS(rf). (5) 

k 

Then, the difference between I(p AB ) and C(p AB ) gives the amount of quantum discord in the system 

D( P AB ) = I{ P AB ) - C(p AB ) = S(p A ) + S min - S(p AB ). (6) 

The main difficulty, for several qubits as well as qudits systems, lies in performing the minimization of 
the conditional entropy. This explains why there is no straightforward algorithm to compute explicitly 
quantum discord for mixed states. Only partial results are available. They were obtained for some 
special forms of the so-called X-states [281 ED Ell EH] . 

3 Amplitude damping 

The beam splitter offers a simple way to probe the quantum nature of electromagnetic field through 
simple experiments. The study of entangled states has revived interest in this device. Many authors 
have considered the behavior of quantum states when passed through a beam splitter [361 139], Recently, 
a quantum network of beam splitters was used to create multi-particle entangled states of continuous 
variables [40] and also multi-particle entangled coherent states [JI]. It also provides, as mentioned 
in the introduction, a simple way to model the amplitude damping related to the absorption of 
transmitted photons in a noisy channel. 

3.1 Fock state inputs 

The beam splitter is an optical device with two input and two output ports that, in some sense, 
governs the interaction of two harmonic oscillators. The input and output boson operators are related 
by a unitary transformation which is an element of the SU (2) group defined by 



B{6) = exp 



— [a 1 a 2 — a{ a 2 



2 



(7) 



The objects and a ; (I = 1,2) are the usual harmonic oscillator ladder operators. The reflection 
and transmission coefficients 



t = cos — , r = sm- 8) 
2 2 w 

are defined in terms of the angle 9. The operator B is actually acting on the states |ni,?i2) °f the 
usual two dimensional harmonic oscillator. If the input state is \ni,n2), then the B action leads to 
the following Fock states superposition 

B\ni,n 2 ) = y~] (rai,m 2 |g|rai,n2)|mi,m 2 ) = ^ B™\i£ 2 \mi,m 2 ) (9) 

mi,m2 mi,m,2 

and in general the output is a two-particle entangled state. On the other hand, the action of the 
unitary operator B on the state |n, 0) gives 



B\n,0) = (l + |g|T"]Pr 7= V = \n-m,m) (10) 
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where the new variable £ = r/t is defined as the ratio of the reflection and transmission coefficients 
of the beam splitter under consideration. Then, the output state (jlOp turns out to be the SU(2) 
coherent state associated with the unitary representation labeled by the integer n. This method can 
be extended to a chain of k beam splitters to generate SU(k + 1) coherent states labeled by the 
reflection and transmission parameters (see for instance [42 \ |4"3"] where similar notations were used) . 
It is important to stress that the generation of coherent states using beam splitters requires input 
radiation state with fixed number of photons. The experimental production of such interesting and 
highly non classical states has been investigated during the last decade (see [H] and references therein). 
Recently an important experimental advance was reported by Hofheinz et al in |45j. They gave the 
first experimental demonstration for generating photon number Fock states containing up to n = 6 
photons in a super-conducting quantum circuit. 

3.2 coherent state inputs 

To describe the photon loss, we usually assume that some of the coherent field is lost in transit via a 
beam splitter. The coherent states enters one port of the beam splitter and the vacuum, representing 
the environment, enters the second port. After transmission some information encoded in the coherent 
states is transferred and the remaining amount of information is lost to the noisy channel. To find 
the final state after transmission, we should first evaluate the action of the operator B{6) on the 
transmitted state. Here we shall be interested in the superpositions of the form 

|Q Q > = _L=(a| -«> + %*)) (11) 

V N i a ) 

where \a\ 2 + \b\ 2 = 1 and N(a) is a normalization factor given by 

N(a) = l + e- 2 ^\ab* +a*b). 

So, we start by evaluating the beam splitter action on the bipartite state \a,0) where the Glauber 
coherent-state \a) is expressed in the Fock (number) basis as 



oo 

a' 

\a) = e 2 

n=0 



This yields 

B(0)|a,O} = \at,ar). (13) 
It follows that the action of the operator B{9) on the state \Q a ) gives 

\Q)t = . (a\ - at, -ar) + b\at, ar)), (14) 
V N ( a ) 

a state describing the quantum field and loss modes (environment). The final state is then obtained 
by performing a partial trace over the loss mode I. We get 

oo 

P = Y J MQ)tt{Q\n) i . (15) 

n=0 
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A straightforward calculation gives 



1 + c 



1 



(a I — at) + 6|ai)) I x h.c 



1 -c 



: (a| — at) — b\at)) 



x h.c. 



(16) 



where h.c. stands for Hermitian conjugation and c is the the coherent states overlapping c = e 2r2 M 2 . 
The last equation can be also written as 



N(at) 
N(a) 



\{l + c)\Qat){Qat\ + \{l ~ c)Z\Q at )(Q at \Z 



(17) 



where Z is the Pauli Z-operator defined by 

1 



Z\Q 



at 



Z 



(a\at) + b\at)) 



VN(cd) 



(a\at) — b\at)) 



.VN{cd) 

which produces a phase flip in the qubit basis. This implies that the transmission of quantum infor- 
mation encoded in optical coherent states suffers from two main types of error: the reduction of the 
amplitude of the coherent state and the phase flip generated by the application of the Pauli Z-operator. 

3.3 Two mode coherent state inputs 

The above considerations can be extended to two mode coherent states of the form 



\Xa,a) = —^==={\/u \a,a) + e'Vl - u> \ - a, -a)) 



(18) 



where 



Ml = 1 + 2^(1 -to) cos 0e" 4 l< 

The action of a beam splitter on the state \xa,a) <8> |0) gives 

1 



\x)t = \Xa,*)®\0) 



{y/to \a, at, ar) + e* \/\ 



a, —at, —ar)). 



The trace over the lost modes gives the density 

P = Tr* \x)t t(x\- 

Here again, one can see that the density p takes the following compact form 



at 



„(1 + c )\Xa,at) (Xa,at\ + ^(1 ~ c)Z\x a ,at) (Xa,at\ Z 



(19) 



(20) 



(211 



The final state is mixed and the amplitude of the second mode is reduced. 



4 Evolution of Bell cat-states correlations under amplitude damping 

The Bell states are very interesting in quantum optics and have been used in the field of quantum 
teleportation and many others quantum computing operations. The experimental generation can be 
realized by sending a cat states of the form | v2o) + | — v2a) and the vacuum into the two input ports 
of a 50/50 beam splitter B(ir/4) (see equation (Tl3|) ) to get 

1 



\B n 



: (|a, a) + I — a, —a)) 



(22) 
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where the normalization factor is defined by 

iV Q = 2(l + e- 4 H 2 ). 

Obviously, the problem of generating Bell states is deeply dependent on the availability of a source of 
cat states. For instance, they can be produced by sending a coherent state into a nonlinear medium 
exhibiting the Kerr effect They can be also generated using a squeezing interaction, linear optical 
devices and photon counters |47^ 148] , Recently, a promising new method to produce cat states was 
proposed by Lund et al. in [39] . The production of cat states especially ones of high amplitude or mean 
number of photons remains an experimental challenge. Considering the fast technical progress and 
the increasing number of groups working in this field, we expect that the generation of cat states (and 
Bell states) is a goal that is achievable in the near future. The available experimental results in the 
literature, obtained with the present day technology, are encouraging. Superpositions of weak coherent 
states with opposite phase, resembling to a small "Schrodinger's cat" state (or "Schrbdinger's kitten"), 
were produced by photon subtraction from squeezed vacuum [50]. Also, the experimental generation 
of arbitrarily large squeezed cat states, using homodyne detection and photon number states (two 
photons) as resources was reported in [50]. Very recently, creation of coherent state superpositions, 
by subtracting up to three photons from a pulse of squeezed vacuum light, is reported in [51]. The 
mean photon number of such coherent states produced by three-photon subtraction is of 2:75. 



4.1 Bell cat-states under amplitude damping. 



Using the results of the previous section, it is simple to verify that under the action of a beam splitter 
on the state (j22j) . the resultant density is 



AB = ^a± 



i(l + c)\B atat )(B aM \ + i(l - c)Z\B aM )(B aM \Z 



(23) 



To study the quantum correlations in this state, a qubit mapping is required. It can be defined as 
follows. For the first mode A, we introduce a two dimensional basis spanned by the vectors \u a ) and 
\v a ) defined by 

\a) = a a \u a ) +b a \v a ) \ - a) = a a \u a ) - b a \v a ) (24) 

where 

|<2a| 2 + \b a \ 2 = 1 |oa| 2 — \b a \ 2 = (—a\a). 
To simplify our purpose, we take a a and b a reals: 

yi +p , yi -v / i \ -2M 2 

a a = -j= — b a = -= — with p = \—a\a) = e 1 1 . 

y2 y2 

Similarly, for the second mode B, a two-dimensional basis generated by the vectors \u a t) and \v a t) is 
defined as 

\at) = a at \u at ) + b at \v at ) | - at) = a at \u at ) - b at \v at ) (25) 



where 



O-at 



y/l+p 1 ' 2 , y / l-p t 

b at = 



V2 



V2 
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The density matrix p AB (|23|) can be cast in the following matrix form 

/ (1 + c)a 2 a a 2 at (1 + c)a a a at b a b at \ 

(1 - c)a? a b 2 at (1 - c)a a a at b a b at 

(1 - c)a a a at b a b at (1 - c)b 2 a a 2 at 

\ (1 + c)a a a at b a b at (1 + c)6 2 6 2 t y 

in the representation spanned by two-qubit product states 

\l) = \u a ) A ® \u a t) B \2) = \u a ) A ® \v at ) B \3) = \v a ) A ® Wat) B \4) = \v a ) A ® \Vat) B ■ 



(26) 



The visual form of the obtained density (|26p resembles the letter X and it is a special kind of the 
so-called called X-states which have been extensively discussed in the literature [52j[53]. The density 
p AB also rewrites, in the Bloch representation, as 

3 



P 



AB 



-(I <g> I + i? 30 CT 3 <8) I + #03 I ® 0-3 + Rii O-i ® <Ji 



where the correlation matrix elements are given by 



R 



03 



p + p 



2-r 



l+p 2 



R. 



30 



i=l 



2p 
l+p 2 ' 



i? 



ii 



V(i-p 2 )(i-p 2 * 2 ) 



-R22 = -p 



l+p 2 

4.2 Quantum mutual information 



l+p 2 



i?: 



33 



l+p 2 



(27) 

(28) 
(29) 



The density p" 4 ' 8 is a two qubit state of rank two. The corresponding non vanishing eigenvalues are 
given by 



A 1 



(l+p r2 )(l+/ +1 ) 



(l-p r2 )(l-p t2+1 ) 
2 + 2p 2 



2 + 2p 2 

It follows that the joint entropy is 

S{p AB ) = -Ai log 2 Ai - A 2 log 2 A 2 . 

The quantum mutual information is given by 

I{p AB ) = S(p A ) + S(p B ) + ]T A, log 2 A,- 

3=1,2 

where p A and p B are the marginal states of p AB , and 

S(p A ) = -X A log 2 A^ - A^ log 2 X A S(p B ) = -X B log 2 X B - X B log 2 X B 

with 

X A_ (1±P? x? _ (1 ±/)(l± P r2+1 ) 



(30) 
(31) 
(32) 

(33) 



2 + 2p 2 



2 + 2p 2 



Reporting (|33|) into (|32|) . the quantum mutual information reads 



- Hi {1+P) 'l ) +H ( {1+pt2){1 +pr2+V 



2 + 2p 2 J V 2 + 2p 2 
where H(x) = —x log 2 x — (1 — x) log 2 (l — x). 



H 



(i+p r2 )(i + p t2+1 ; 

2 + 2p 2 



(34) 
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4.3 Conditional entropy 



After computing the quantum mutual information, we need next to compute the classical correlation 
C(p AB ) defined by @. We consider projective measurements for subsystem A. We follow the proce- 
dure developed in [53]. We remind that the generalized positive operator valued measurement is not 
required. Indeed, as mentioned above, the optimal measurement for the conditional entropy is ensured 
by projective operator [27] . The general form of the SU(2) unitary operator, occurring in ([3]), is 

U = exp(rja + — fja_) exp(i(j)a^) (35) 

where rj £ C and (f) G R. This parametrization allows us to express the quantities defined by 

{<Ti) k = {k\U^<TiU\k), i = 1,2,3 and k = 0, 1 



as 



k 1 — zz _ k z + z 



1 + zz 1 + zz 1 + zz 

where z = — i-^U tan \pf\r\. They can be also written as 

(o"3>fc = (-) fc cos#, (cri) fe = (-) fe sin6»cosv9, <cr 2 ) fe = {~) k sin 9 sin p 

where |e i</3 = —ir). These mean values combined with the expressions of the correlation matrix 
elements (|28p and (j29|) determine the explicit form of the conditional densities ([2]) and the conditional 
entropy ([5]) in terms of the angular variables 9 and p. Indeed, combining the equations d21) , (J3j) , (j26|) 
and (|35|) . it is simply seen that the density operators p B take the following form 



? B = J_( (1 + ^03) + (#30 + #33)(^ 3 )fc Rll{<Tl)k - «#22 Wfc 

" k Pk I Rll{°l)k + *fl22 W* (1 " ^03) + (#30 - «33)(cT 3 )fc 



where 



Pfc =^(1 + ^30 (a 3 ) k ). 



It follows that the conditional entropy given by ([5]) rewrites also as 

S = S(9,<p) = Y,Pk tf(V^l-4de t/9 f) (37) 
fc=0,l ^ ' 

and can explicitly expressed as a function of # and p. Then, the minimization of S can performed 
over the polar and azimuthal angles. Nevertheless, there exists another elegant way to optimize the 
conditional entropy. It is based on the Koashi- Winter relation [26] (see also [19]) as we shall explain 
in what follows. 



4.4 Minimization of conditional entropy 

The Koashi- Winter relation establishes the connection between the classical correlation of a bipartite 
state p AB and the entanglement of formation of its complement p BC . This connection requires the 
purification of the state p AB . In this respect, as p AB is a two-qubit state of rank two, it decomposes 
as 

p AB = A 1 |Vi)(^i| + A 2 |^ 2 )(^ 2 | (38) 
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where the eigenvalues Ai and A2 are given by (|3"Uj) and the eigenstates \tpi) and \tp2) are 

= 1^2 U9U2 ( a a a at\ u a,U at ) + b a b at \v a , V at )) 
V a a a at + K 13 at 

1^2) = —^=====(a a b at \u a ,v at ) + b a a at \v a ,v at )). (39) 

Attaching a qubit C to the bipartite system AB, we write the purification of p AB as 

/ AT|'0i> ® K> + V^M) ® K) (40) 

such that the whole system ABC is described by the pure density state p ABC = \ip){Tp\ from which 
one has the bipartite densities p AB = Trcp ABC and p BC = TtaP ABC ■ Suppose now that a von 
Neumann measurement {Mo, M\} is performed on qubit A (here also we need positive operator valued 
measurement of rank one that is proportional the one dimensional projector). From the viewpoint of 
the whole system in the pure state the measurement gives rise to an ensemble for p that we 
denote by 

s BC = { Pk M BC )} 



where 

Pk 



(Af fc <g>I®I)|^)M 



On the other hand, from the viewpoint of the state p AB , the von Neuman measurement on A gives rise 
to the ensemble for p B defined previously as £ B = {p B ,p B }- It is simple to check that the ensemble 
£ B can be induced from £ BC by tracing out the qubit C, namely 



Pk = Trc 



^ C H<P BC \ 



We denote by E(\4> BC )) the measure of entanglement for pure states. It is given by the von Neumann 



J k 

entropy of the reduced subsystem p B = Trc(|^ c ')(^ c '|) 

E(\<p BC )) = S(p B ). 

It follows that the average of entanglement of formation over the ensemble £ 

E BC = £ Pk E^ BC )) 



-BC _ 

fc=0,l 



coincides with the conditional entropy ©.At this level, it is important to notice that Koachi and 
titer 

of p BC . Consequently, 



—BC 

Winter [26] have pointed out that the minimum value of E is exactly the entanglement of formation 



S min = E(p^). (41) 

This relation simplifies drastically the minimization process of the conditional entropy. Therefore, 
employing the prescription presented in [55] to get the concurrence of the density p BC , one obtains 



E(p BC ) = h(^+ l -^l-\C{pBC)\2^ (42) 
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with 



c{pBC) = y^(l-p^)(l- p 2t- 



(l+p 2 ) 

It is remarkable that the minimal value of the conditional entropy given by (|42p is reached for 9 = tt/2 
and tp = 0. Indeed, using the explicit form of the conditional entropy (|37p in term of the polar and 
azimuthal angles 8 and (j>, one can verify that 

s min = s(e = ^,4> = o). 

It is important to stress that the special values 9 = and <j) = n/2 coincide with ones obtained in |18| 
(see also [SB]) for optimal measurements to access quantum discord of two-qubit states. According to 
the equation (J3|), the classical correlation is 

<y ») = h( 1 + + C )- b( l - + l - £± p2 + P "T + p2 " + 1 (43) 
^ 7 V 2 2 l+p 2 / V 2 2 l+p 2 / 

and using the definition ([6]), the explicit expression of quantum discord reads 



VF ; V 2 1+P V V 2 2 X +P / V 2 2 !+P / 

Note that for r = 0, the density state p AB ([23]) reduces to the pure density Bell cat-states (|22|) and 
the the quantum discord (|44j) gives 

'1 + _P. 

which is, as expected, exactly the entanglement of formation E(\B a ^ a ) (B a ^ a \) . Indeed, it is simple to 
check that the concurrence of the pure Bell cat-state \B a ^ a ) is 

1-p 2 
l+p 2 

which coincides for r = with the concurrence of the state (|23|) given by 



D(\B a>a )(B a , a \) =H(- + ^3 ) (45) 



C(\B a>a )) = (46) 



C(P AB ) = " Wy W:; - ' ■ (47) 



/V(l -p 2 )(l ^) 
l+p 2 



where we used the definition of the concurrence for mixed states introduced in [55]. It is also simply 
verified from (|44p that for r = 1, the quantum discord vanishes and the concurrence (|47|) is also zero. 

The behavior of quantum discord (|44p versus the overlapping p and the reflection coefficient r is 
plotted in the figure 1. As seen from this figure, for p fixed the quantum discord decreases as we 
increase the reflection parameter r of the beam splitter and the maximum is obtained in the limiting 
case r = 0. This indicates that the noisy channel, inducing the decoherence effects, renders the system 
less correlated and subsequently the quantum discord decreases. In other hand, for a fixed value of 
r, the quantum discord starts increasing, reaches a maximal value and decreases after. To see this 
behavior, we plot in the figure 2, the quantum discord versus the overlapping for Bell cat-states for 
different values of the reflection coefficient r. In particular, under the action of a 50:50 beam splitter, 
the maximal value of quantum discord is reached for p ~ 0.4. Clearly, for r = 0, the maximum of the 
quantum discord or entanglement of formation is attained for p = 1. From the figure 2, it is easily seen 
that the value of the overlap p, for which the quantum discord D is maximal, increases for increasing 
values of r. 
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FIG. 1: The pairwise quantum discord D versus the overlapping p of Glauber states and the reflection 

parameter r 2 . 



D 




P 



FIG. 2: The pairwise quantum discord D versus the overlapping p for different values of r 

5 Evolution of quantum correlations under dephasing channel 

In this section, we investigate the dynamics of bipartite quantum correlations (entanglement and 
quantum discord) of damped Bell cat-states, given by the density p AB (TZBj) . under dephasing channel. 
To describe conveniently the effect of this channel, we use the the Kraus operator approach (see for 
instance [5])- In this approach, the time evolution of the bipartite density p AB = p AB (0) (|26p can be 
written compactly as 
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where the so-called Kraus operators 



E„Jt) = EJt)®E v {t) 



i. 



The operators E^ describe the one-qubit quantum channel effects. For a dephasing channel, the 
non-zero Kraus operators are given by 



E = diag(l, y/T^) E 1 = diag(l, ^7) 
with 7 = 1 — e _rt and V denoting the decay rate. This gives 



/ 



„AB 



(t) = ~ 







c)a 2 a b 2 at 



(1 + c)a 2 a a\ 



" ri (l + c)a a a at b a b c 
The entanglement in this state is measured by the concurrence 







(1 

-rt(i. 



-rt. 



1 - c)a a a at b a b a t 



V 



J at 



c)a a a at b a b at 




:i 



c)bW 



at 



e~ rt (l + c)a a a at b a b at \ 



(l + c)b 2 a b 2 at 



C = 2 max{0, Ai(t),A 2 (t)} 



where 



Ai(i) = —a a a a tb a b a t 

1\ rv 



(l- 7 )(l + c)-(l-c) 
Since A 2 (t) is non positive, the concurrence is 



A2(i) = —a a a a tb a bat 



N, 



(l- 7 )(l-c)-(l + c) 



cm _i V(i-p 2 )(i- 



2t 2 ^ 



e - ri (l+/')-(l-p^) 



(49) 



for 



i < t 



i[ ln(l+p r 



■ln(l-P r )]■ 

In this case, the system is entangled. However, for t > to, the concurrence is zero and the entanglement 
disappears., i.e. the system is separable. This shows clearly that under dephasing channel, the 
entanglement suddenly vanishes. Note that the bipartite system under consideration is entangled 
in the absence of external noise. Indeed, for t = 0, the concurrence (|49p reproduces (|47p which is 
non zero except in the limiting case p — > ( \a\ — > +00) or r = 1. The phenomenon of total 
loss of entanglement, termed in the literature "entanglement sudden death" [57], was experimentally 
confirmed under some specific conditions [58]. As the concurrence vanishes after a finite time to, it is 
interesting to ask what happens to quantum discord. The explicit expression of the amount quantum 
discord in the state p AB (t) can be obtained following the method described in the previous section. 
But, here we need only to know if the state p AB {t) has vanishing quantum discord. As the density 
p (t) has also the form of the letter X, one can use the criteria, classifying the so-called X states with 
vanishing quantum discord, discussed in [59]. Therefore, according to this criteria, p AB (t) has zero 
quantum discord if and only if p — > 0. The quantum discord in the state p AB (t) is in general nonzero 
even when entanglement suddenly disappears. This gives a special instance of separable quantum 
states for which the quantum discord is non zero. This agrees with the commonly accepted fact that 
the quantum discord is a kind of quantum correlations which goes beyond entanglement and almost 
all quantum states have non vanishing quantum correlations 
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6 Concluding remarks 

To close this paper, let us briefly summarize the main results. Our effort was devoted to investigate 
the evolution of quantum discord of Bell cat-states under amplitude damping channel. We have 
derived the analytical expressions for the classical correlation and quantum discord. We discussed the 
usefulness of the Koashi- Winter relation to determine a closed form of quantum discord present in the 
system. We also discussed the evolution of the transmitted Bell cat-states under a dephasing channel. 
This provides a special instance for separable mixed quantum states with non vanishing quantum 
discord. The method used in this paper, to derive quantum discord, constitutes an alternative way to 
compute analytically classical correlations and quantum discord. It can be applied easily in evaluating 
the quantum discord present in the other Bell cat-states: 

\Ba,a) ~ |«, a) - | - a, —a) 

Wt-a) ~ l a ; ~ a ) + I - «, «) 
\B~a-a) ~ \ a i ~ a ) ~ I ~~ a i a ) 

after passing trough a amplitude damping channel. In addition, the analysis presented here for the 
evaluation of quantum discord is readily extended to more general systems, including squeezed states, 
SU(2) and SU (1,1) coherent states and so on. Finally, it will be interesting to compare the quan- 
tum discord in bipartite coherent states with its geometrized version usually called in the literature 
geometric quantum discord [ST]. Further thought in this direction might be worthwhile. 
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